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We apply linear perturbation theory to the study of the universality and criticality first observed 
by Choptuik in gravitational collapse. Since these are essentially nonlinear phenomena our attempt 
is only a rough approximation. In spite of this, universal behavior of the final black hole mass is 
observed with an exponent of 1/2, slightly higher than the observed value of 0.367. The universal 
behavior is rooted in the universal form that in-falling perturbations on black holes have at the 
horizon. 



I- If one considers the collapse of a spherically symmetric scalar field in general relativity, two possible end-results are 
expected. Either the initial data manages in the time evolution to create a big enough mass concentration such that 
a black hole forms or whatever scalar field was present is dispersed to infinity, leaving a flat spacetime as a result. 
If one considers one-parameter families of initial data, there will be in the parameter range two distinct regions, one 
\ corresponding to the formation of a black hole and another to dispersion to infinity. Christodoulou and others 

had raised the question of what was the behavior of the final black hole mass as a function of one such parameter. 
For values of the parameter below a "critical" value for which black holes start forming, the final mass is zero, since 
no black hole forms. When the critical value is reached, is there a "jump" in the mass, i.e., is there a minimal mass 
for which one can form a black hole? This was highly unlikely, since the problem has no scale. Choptuik Q has 
confirmed this through numerical simulations. He finds no "mass gap" , but actually finds the final black hole mass is 
&\ • a universal function of the parameter value of the initial data, 

o : 

^ . Mbh = k(P ~ P*) 13 (1) 

G\ ■ 

where k is a constant that may be different for different one-parameter families of data, p is the parameter value 
chosen for the initial data and is the critical parameter value at which black holes start to form. The exponent is 
universal (i.e. family-independent) and was experimentally determined by Choptuik to be f3 = 0.367. This value was 
bX). later confirmed with a different (less sophisticated) kind of code Q . 

In a separate set of experiments, Abrahams and Evans [^) considered the implosion of an axially symmetric grav- 
• i-H , itational wave (in vacuum, no scalar field present). They find a similar behavior with the same exponent. Finally, 
Coleman and Evans || considered the collapse of a radiative fluid p = p/3 and also found a similar behavior with the 
5h I same exponent. 

Little theoretical progress has been made towards explaining this universal behavior and in particular the value of 



the exponent. Several authors have noticed that there exist in the literature exact self-similar solutions for the 
collapse of scalar fields. These solutions are not quite geared towards describing the effects observed by Choptuik, 
since they are not asymptotically flat. In spite of this they yield a critical behavior with an exponent of 1/2. An 
exponent of f/2 has been also observed in a rather unrelated context, the collapse of f + f dimensional black holes 
with quantum corrections by Strominger and Thorlacius flO|| . 

The purpose of this note is to explore the application of perturbation theory to these phenomena. Our approach 
is crude, and the results can only be taken as suggestive. It is remarkable however, that we find a universal exponent 
that is close to the observed value (1/2 instead of 0.367), and that the root for the universality is the same as that of 
the no-hair behavior of black holes: the universal form that fields falling into the horizon take (quasinormal ringing) . 

Suppose one prepares a set of initial data for the Choptuik experiment, with a certain value of the parameter p\ 
(although it is not quite needed, the reader can imagine a Gaussian of amplitude p\ for concreteness) . Let us also 
assume p\ > p* and therefore, according to Choptuik, a black hole forms, with mass, 

A/i = kfa - p.) 13 . (2) 

Suppose now a second set of data is prepared, of parameter value pi > pi, with p2 differing little from pi, i.e., 
Pi — Pi < Pi — P*- Again, we form a black hole of mass, 
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M 2 = k(p 2 -p*f 



(3) 



and we are assuming the two sets of data are in the same family, so the constant k is the same. We now would like 
to view the second set of data as a "perturbation" of the first set. Assuming their parameter separation to be small, 
i.e. Ap = p2 — pi « Pi — p*, we can write for the final mass difference AM = M 2 — Mi, 

AM = kf3( P2 - Pl )( Pl - p*) 13 - 1 (4) 

which can be written as, 

i 3-1 

AM = k* pApMpf . (5) 

Therefore, if there would be a way of estimating the excess mass that falls into a black hole due to an increase of 
the parameter value of the initial data, one could estimate the exponent (3. Alternatively, if one could estimate AM 
as a function of the black hole mass Mi one could also get a value for the parameter (3. 

We will now provide such an estimate using perturbation theory. We start by making a rough assumption: the 
spacetime geometry for the collapse of the first initial data set can be approximated by the spacetime geometry of 
a static, ever-existing, Schwarzschild black hole of mass M\. This is an uncontrolled and vague approximation. One 
expects it in general grounds to be good if p\ is far away from criticality, so one forms a black hole with most of the 
initial data mass and there is little outgoing radiation. Also one should realize that approximating a time-dependent 
spacetime by a static one can only work in the region where the time-dependent spacetime is not rapidly changing. 
In fact, numerical studies of perturbations on dynamical backgrounds have shown that the approximation can work 
remarkably well 0]. However, there is no real justification to apply this approximation in the case we are interested 
since the attention should be focused near criticality, when the background differs considerably from an ever-existing 
black hole. What follows therefore should only be taken as suggestive and speculative. 

If one assumes the spacetime determined by the initial data set with parameter p\ is a Schwarzschild black hole 
of mass Mi, estimating AM becomes a problem in perturbations of black holes: given a perturbation incident on 
a black hole, how much mass falls into the hole and how much mass is radiated away? To estimate this we note 
that perturbations on black hole backgrounds are determined by the Zcrilli function \I/ (to summarize, one takes 
9nu = Schwarzschild^ + and finds that all the information in h^ v can be parameterized by a function \E', see pd| ] 
for details). This function satisfies a Klein-Gordon- like equation coupled to a potential called the Zerilli equation. 
The "radial" variable in this equation is a "tortoise" coordinate r„ in terms of which the horizon is at r„ = — oo and 
i Q is at r„ = oo. The question is therefore, given a perturbation incoming from r* — oo, how much energy makes it 
to = — oo? Scattering in the Zerilli potential has been well studied, and one knows that for a generic incoming 
perturbation from = oo one gets a damped oscillation known as quasinormal ringing at the horizon (r* — — oo), 

-^Horizon = A exp(-w;i) sin(w r t) (6) 

where the quasinormal frequencies scale as inverse powers of the background spacetime mass, 

constant 

w - = ~mT (7) 

and the coefficient A depends on particular details of the initial perturbation. A study of how different perturbations 
generate different values of A is found in Price and Sun . 

The total energy falling into the black hole is proportional to the integral J dt^ olizon . The dependence of this 
integral on Mi is determined completely by the behavior of the quasinormal frequencies and gives as a result, 

AM - Mf 1 (8) 

which implies (3 = 1/2. 

Let me end this note with some remarks. 

• The above explanation, although rough, provides a link between universality and no-hair behavior in gravitational 
collapse. 

• In spite of the coarse approximation, the exponent comes out with a value close to the experimentally measured 
one. It seems that the fact that the value is the same as in the exact solutions or the Strominger-Thorlacius |lC| ] 
model is a mere coincidence, since none of the features of those examples has been incorporated in this discussion. 

• The above mechanism predicts the same value for the exponent for the scalar field studied by Choptuik, the 
gravitational waves studied by Abrahams and Evans and the fluid case studied by Coleman and Evans, since it only 
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depends on the relation of the quasinormal frequencies to the background spacetime mass, which does not depend on 
the details of the perturbations. 

The outstanding question is: is there any way to refine the approximation in order to give the correct exponent? 
The following remarks are in order. 

• The dependence of quasinormal frequencies on the background mass is a robust feature, that seems unlikely to 
change if one modifies the scenario. It is possible that a dynamical background could alter the picture of quasinormal 
ringing at the horizon and therefore alter the exponent. This would bring into question up to what extent is the 
link between universality and baldness a mere artifact of our rough approximation. If quasinormal ringing were still 
present in a dynamical background, in order to predict the right exponent, the frequencies would have to depend 
non-analytically on the background mass. This seems quite unlikely, even on dimensional grounds. Traschen |p_3[| has 
found some non-analytic behavior in perturbations of Reissner Nordstrom. Perhaps there is a connection. It could 
also happen that in the non-static case there is a residual dependence of the coefficient A on Mi , which in the static 
case does not appear. This could allow non-analytic behavior of the energy without a counter-intuitively non-analytic 
behavior of the quasinormal frequencies. 

• Insight into perturbations of dynamical backgrounds could be gained by studying perturbations of the Vaidya 
metric. This would provide a dynamical background and could help gain intuition on some of the issues raised in the 
previous point. 

• Finally, in the case of the collapse of radiative fluid, studied by Coleman and Evans, the exact critical solution 
is known analytically. Studying perturbations of it should certainly allow insight into the critical exponent. Evans is 
currently pursuing this point, with a rather different approach than the one presented here. 
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